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Wideband pulse propagation: single-field and multi-field approaches to Raman interactions 
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We model the process of ultra broadband light generation in which a pair of laser pulses separated by the 
Raman frequency drive a Raman transition. In contrast to the usual approach using separate field envelopes for 
the different frequency components, we treat the field as a single entity. This requires the inclusion of few-cycle 
corrections to the pulse propagation. Our single-field model makes fewer approximations and is mathematically 
(and hence computationally) simpler, although it does require greater computational resources to implement. 
The single-field theory reduces to the traditional multi-field one using appropriate approximations. 
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I. INTRODUCTION 

An important aim of current wideband Raman experiments 
is to efficiently generate few-cycle pulses ^ H,^]. If driven 
strongly enough, the two-photon Raman transition modulates 
the incoming field by adding sidebands separated by the tran- 
sition frequency. Wideband fields are generated as these side- 
bands generate sidebands of their own (and so on); a wide 
comb of frequency components separated by the transition 
frequency is generated in this way. If a scheme can be imple- 
mented that adjusts the phases of each component appropri- 
ately, then few- or single- cycle optical pulses can be obtained 
(see e.g. ||]). 

In standard theoretical treatments of the Raman process, the 
field is split into frequency components centred on the teeth 
of the frequency comb. This approach has the advantage that 
the components can be modeled reasonably well with slowly 
varying envelopes, but of course it has the disadvantage that 
one needs to keep track of a large number of components. In 
this paper, we study an alternative approach in which the field 
is treated as a single entity rather than being split into pieces. 
Note that this approach is distinct from methods based on the 
direct solution of Maxwell's equations such as FDTD (finite 
difference time domain) [H] or PSSD (pseudospectral spatial 
domain) [|[|. Our single-field is based on a second-order wave 
equation, and uses a convenient choice of carrier to define a 
field envelope. As we will demonstrate, the latter technique 
offers significant advantages over the traditional multi-field 
formalism. 

To provide a context for the discussion, we consider exper- 
iments such as those of Sali et.al. 0), where the Raman 
transition is driven near-resonantly by a pair of intense pump 
pulses about lOOfs long; compared to the transition frequency 
of about 130THz, the spectra of each pump pulse (and hence 
the generated sidebands) are relatively narrow. This means 
that a multi-component model is still not unreasonable, even 
if numerical considerations might demand that the aiTays used 
to store these spectra overlap in frequency space. However, if 
we were to move to shorter pump pulses, or to a single (much 



shorter) pump pulse with enough bandwidth to efficiently ex- 
cite the transition, we would reach the regime where the teeth 
of the frequency comb significantly overlap. At this point, 
one would be forced not only to move from a solution based 
on the Slowly- Varying Envelope Approximation (SVEA) to a 
more accurate model such as the Generalized Few-Cycle En- 
velope Approximation (GFEA) but the utility of mul- 
tiple field components would in any case become question- 
able. This provides the motivation for the present work, since 
our model can extend the regime in which the advantages of 
envelope-based methods can be utilised; it also turns out to 
be more versatile, placing fewer restrictions on the kinds of 
Raman media that can be easily described. 

In this paper, we will construct a single-field model which, 
in most other respects, closely parallels the approach to wide- 
band Raman generation adopted by Hickman et. al. JlO|]. 
A key feature of the single-field model is that the coupling 
constants oscillate at the Raman frequency, and it is this that 
impresses the sideband modulation on the propagating field. 
Since the field is now not only wide-band but contains signifi- 
cant sideband components (i.e. distinct sub-peaks, as opposed 
to a broad featureless background), the field envelope is no 
longer slowly- varying and must therefore be propagated using 
the GFEA. This necessity can be demonstrated by comparing 
the results of the single-field model with those of a multi-field 
counterpart. 

The paper is organized as follows: section || outlines the 
derivation of the single-field Raman theory, section |l^ shows 
how to reduce it to a standard multi-field version, and section 



IV 



applies the theory to practical situations. In section ^ we 
discuss some of the issues relating to our Raman model and 
its numerical implementation, and finally section ^ contains 
our conclusions. 



II. SINGLE-FIELD RAMAN THEORY 

We start by considering the wave function \|/ of a single 
molecule (e.g. H2) and the electric field E, and write the 
time-dependent wave function by expanding it in terms of the 
eigenfunctions in the field-free (i.e. £ = 0) case. This means 
we can get the expansion coefficients by solving for an ef- 
fective Schrodinger equation that contains a two-photon Rabi 
frequency created by means of an interaction term based on 
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a field-dependent dipole moment. We assume a dispersion- 
less medium and write all equations in terms of position z and 
retarded times t = tiab — z/c. Here we follow the method of 
Hickman, Paisner, and Bischel Jiofa (HPB), but we use only a 
single E field rather than multiple components. Note that HPB 
use Gaussian units, so there may appear to be inconsistencies 
when comparing our formulae (in S.I.) to theirs. 

We denote the known molecular eigenfunctions of the un- 
perturbed Hamiltonian Hq as \n), and their corresponding en- 
ergies hW„. We want to obtain the solution to 
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where d is the electronic dipole moment operator and the c„ 
are a set of complex probability amplitudes. 

We now replace the electric field E with a carrier-envelope 
description, but, unlike HPB, we use only a single component 
centred at a frequency of coo, rather than a set indexed by an 
integer j. The envelope and carrier for the field is: 

E = Ae-'f^o'-'^o'^+cc, (4) 

and, following the standard procedure of assuming the co- 
efficients c, are slowly varying, discarding terms at multiples 
of the carrier frequency, and simplifying, we eventually reach 



ih- 



dt 



(5) 



where a„j = -exp[-iW,„f]^<i„,<ii^ 



COS 



(6) 



The a„j coupling parameters oscillate because, in contrast 
to the HPB derivation, there is no frequency difference be- 
tween field components to cancel with the Raman transition 
frequency. We now take the indices 1 and 2 to correspond 
to the two states involved in the Raman transition of inter- 
est; these will be the and 1 vibrational (or perhaps rota- 
tional) levels of the electronic ground state. Indices 3 and 
above will correspond to (quoting HPB) "translational motion 
on higher electronic states". Since we are interested only in 
the Raman transition, we specialize the above equations for 
the coefficients c„, calculating ci and C2 only, and assuming 
that the dn = {l\d\2) dipole moment is zero. This means 
we will only be including transitions between indices 1 and 2 
that go via one of the higher states j > 3, since we still allow 
dij, d2j ^0 VJ > 3. Further, we solve for the coefficients for 
the higher states in terms of ci and C2, in an adiabatic approx- 
imation justified when ci and C2 vary only slowly compared 
to the exponential terms. 

When converting the equations for ci, cj into Bloch equa- 
tions, we make the same approximations as HPB: keeping the 
energy separations for all transitions greater than that of the 
1^2 transition, and ignoring all the higher vibrational (or 
rotational) states. Thus we can write 



at2-0C2i ~ 0, 



(7) 
(8) 



Here co/, is the Raman transition frequency, and 5' is a phase 
factor that ensures that the coupling constant /' is real valued. 
This /' will be used to replace an + ot^i- We also get a Stark 
shift term - 



%' = «ii-«2V 
We define pi2 = cic^ and w — C2C2 — CjCi, so that 
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.(^^li-^2|APp,2 + .^2|A|^w, (10) 



i^2|A|2p,2-;— ^2|A|2p*2. 



(11) 



Finally, we insert decay terms y,, and introduce coj^ = co^, — 
A. This A allows for arbitrary rotations of the polarization, 
P12 = pj2exp(— ;Af — ;5'). Eqns. ( 10, ll]) governing the re- 
sponse of the medium to the applied fields now become 



^tP'n = (-72 
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The parameter A should be chosen to optimise computational 
accuracy by making the dynamics as slowly-varying as possi- 
ble. For example, if the field contained two frequency compo- 
nents that were slightly detuned from the Raman frequency, 
we might use A to compensate for the resultant beating. In 
general, A is most useful in the multi-field model discussed in 
the next section. The complementary part that specifies how 
the field responds to the polarization of the Raman transition. 
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Here the 1 +;3,/coo in eqn.(|l5|) is (with dt = djdt) the 
lowest-order approximation to the GFEA few-cycle propa- 
gation corrections ^, which is equivalent to the SEWA 
(Slowly Evolving Wave Ap proximation) correction derived 
by Brabec and Krausz [pT|]. Although the full form is not 
included for reasons of brevity, it could easily be introduced if 
the extra accuracy was desired; indeed we routinely use it in 
our simulation codes. It is independent of the Raman deriva- 
tion presented here, since it is a field propagation effect. The 
full form of the few-cycle prefactor (and various expansions 
thereof) has already been reported in [||, ^ . 

A detailed derivation of this single-field Raman theory can 
be found in JI^ 

We solve these equations numerically using a split step 
method, where we treat the nonlinearity in the time domain, 
and the dispersion in the frequency domain. To include dis- 
persion in a time domain equation like eqn.(|l5]) requires ei- 
ther additional time derivatives (as in ^) or a convolu- 
tion over a time-response function which is an A^^ operation. 
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However, handling dispersion in the frequency domain is both 
conceptually simpler (since it simply amounts to a frequency- 
dependent phase evolution), and more computationally effi- 
cient because it is an NlogN process. 

The validity of the approximations used in deriving our 
Bloch equations will obviously depend both on the details of 
the chosen Raman medium and/or transition, and on the num- 
ber of Stokes and anti-Stokes sidebands we wish to describe. 
Since in the experiments of [|l], ^ ^ ^ Q the emphasis was 
on a single Raman transition, a simple Bloch model is clearly 
appropriate, and indeed our approximations differ little from 
those of other theoretical approaches (such at that of HPB). 



III. MULTI-FIELD RAMAN THEORY 

The single-field Raman model can be converted into a tra- 
ditional multi-field model as developed in e.g. HPB [ [l0| ] or 
Syed, McDonald and New [|l5]| by replacing the field enve- 
lope with a sum of multiple envelopes using carrier exponen- 
tials spaced at the Raman frequency. When doing this, we 
will only get the correct multi-field form if few-cycle (either 
SEWA or GFEA) corrections to the field evolution part of the 
theory are applied to the effective polarization caused by the 
Raman transition. 

Since the single-field evolution equation (eqn.(p3|)) uses an 
envelope A that is based on a carrier (see eqn.(Q), the single- 
field envelope A is replaced with A ^'s at frequency (aj = COq + 
j'cOi and wavevector kj — k{(Oj). The single-field envelope in 
terms of the new Aj's is 



A = Y,AjSxp[-i{(o'jt~k'jz)] , 



(17) 



where co'- = 03 j — coo, and k'j = kj — k{(0o) ~ kj — ko. 

The equations for and w describing the Raman tran- 
sition result from a simple substitution of eqn.(|l7]) into 
eqns.(p^ pj] ), followed by a rotating wave approximation 
(RWA) to remove non frequency matched terms. They are 

a,p'i2 « (^-j2 + iA + ig'Y.2A*A^p[2 

+ 2if ^ 2AjA*_ 1 .w.e-'^ .e+'i''i-''i-')\ (18) 

dtw — — Yi (w — w,-) , 

+2if (2A*A,+ip;2e'";/ -2A,A*+ip'iV-'"^'') • (19) 

Quite a lot of physics has been removed by the RWA ap- 
proximation, although it is a very reasonable one except in the 
very wideband limit. For example, the effects of next-nearest 
neighbour field components have been ignored, as have all 
more distant field-field interactions. In the next-nearest neigh- 
bour case, the dropped terms would impose a rapid co/, oscil- 
lation onto the polarization pi2, which would in turn tend to 
impose sidebands at ±03^ onto each field component. It is 
reasonable to ignore such sidebands in the narrowband limit 
used for most applications of a multi-field Raman theory; but. 



in principle one might extend a multi-field theory to include 
them by inventing a scheme to apply the sidebands to the field 
component with which they are in nearest resonance. 

Extra factors of 2 have appeared in eqns.([l8|, [l9| ) because 
the multi-field equations start with double summations that 
give pairs of terms that can be reduced to one in the remaining 
single summation. 

Finally, we need to insert the few-cycle correction to the 
polarization term, because the (j ^ 0) sub-envelopes A j have 
an I j(Obt time dependence that cannot be neglected. The po- 
larization correction terms are just the result of applying the 
first-order correction {i/(x)o)dt to the A{t)X{t) from eqn.(15). 
The j-th polarization correction term is then 



ocoap 

{Pi2^./+i exp [+iikj+i - kj)z - iAt\ 



280C0 



-p'i*2A/-i exp [+iikj-i - k'j)z + lAt] } 
i{kj^ko)Aj, (20) 



and differs only from the standard polarization term in that G)'j 
appears in place of coo- The two terms can then be straightfor- 
wardly summed, and since co,- = cOo + CO^, from eqns.(|l5[ 
p^, we get 



dzAj{t) = ^^^^ {p'nAj+i exp - ~ 'Af] 

+Pi*2A;_i exp [+iik'j_i - k'j)z + lAt] } 
-i{kj-ko)Aj, (21) 

where the lAt terms arise because of our rotation of the frame 
of reference of p[2. The residual kj — ^0 terms result from a 
difference in the k frame of reference between the our multi- 
field derivation and the standard one. 



IV. EXAMPLE APPLICATIONS 

We now use the single-field (GFEA) model to simulate an 
experimental situation. Fkst we compare the results to their 
multi-field counterparts, demonstrating the relationships be- 
tween the two methods, and showing them to be in good 
agreement, as expected for the chosen pulse lengths. Second, 
we contrast our model with an (inaccurate) single-field SVEA 
model, in order to highlight the role of the few-cycle propa- 
gation terms. The bulk of the code used was the same for all 
simulations, as it contains options to switch from a single to a 
multi-field case, and to switch GFEA corrections on and off. 

Figure ^ shows a set of results for a pair of pump pulses 
traveling though 9cm of H2. This corresponds to a simula- 
tion of an experiment where the pulses pump the 1st vib(ro) 
level in molecular H2 (at 4155cm^^ i.e. ^ 126THz), as in 
the transient-regime experiments of Sali et.al. [Q, ^]. In these 
experiments, typical pulses might be 70fs and 250fs wide at 
800nm (30/jJ) and 600nm (120/jJ) respectively, and the comb 
of Raman sidebands generated are narrow and well separated. 
A Cauchy-type dispersion curve for H2 is incorporated into 
the simulations. In our simulations, we use the smaller widths 
of 17.5fs and 62.5fs, which broadens the spectral peaks (to 
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about 57THz and 16THz respectively) and makes the stan- 
dard multi-field approach less practical. The figure compares 
three data sets - (a) single-field GFEA simulation, (b) multi- 
field simulation, and lastly (c) single-field SVEA simulation 
(i.e. without any few-cycle propagation corrections). 

There is good agreement in the heights of all the spectral 
peaks between the two exact simulations (single-field GFEA 
fig. 1(a) and multi-field fig[l|(b) ); even the details in the wings 
of the first anti-Stokes peak (at about / ~ 0.5) are replicated. 
Those in the wings of the second anti-Stokes peak (at about 
/ = 0.63) are not well replicated; however, the features in 
question are about three orders of magnitude weaker than the 
peaks, and the two simulations are not equivalent because the 
multi-field theory does not include next-nearest neighbour in- 
teractions. 

The comparison between fig. |l](a,b) and the single-field 
SVEA simulation fig |l](c) is also instructive. Although it does 
reproduce the character of the single-field GFEA spectra in 
many ways, the peak heights do not agree - a fact that is more 
apparent on a linear scale than a logarithmic one. In terms of a 
multi-field model, we can say that without the GFEA correc- 
tions, the prefactor of the polarization term does not pick up 
its correct frequency dependence, so the Stokes lines are arti- 
ficially enhanced, and the anti-Stokes artificially suppressed. 

Figure ^ shows a set of results from a single lOfs probe 
pulse at 397nm, traveling though 9cm of previously polar- 
ized D2. This corresponds to the probe stage of an experiment 
where the gas had been prepared using a pair of nanosecond 
fields resulting in a medium polarization of pi2 = 0.025; on 
the 2993.57cm^' (~ 90THz) vibrational transition, e.g. as in 
the experiments of Gundry et.al. [pj|], who use a longer probe 
pulse of about 150fs. A Cauchy-type dispersion curve is in- 
corporated into the simulations, but in the absence of good 
dispersion data for D2, we use that for H2 as it should be a 
good match. Note that although the polarization initial con- 
dition is fixed, our simulations do incorporate the response of 
the polarization to the probe pulses. The main spectral peaks 
agree well in the multi-field and single-field GFEA simula- 
tions, although as before the results differ at the edges where 
the intensities are very small compared to the main features. 
As for the previous situation, in the single-field SVEA simula- 
tion the Stokes and anti-Stokes lines are artificially enhanced 
or suppressed. 
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FIG. 1: Transient Raman generation using 17.5fs and 62. 5f pump 
pulses as described in the text. Here we compares three simulation 
results: (a) single-field GFEA simulation, (b) multi-field simulation, 
and (c) single-field SVEA simulation. The dashed lines help com- 
pare the relative heights of the first Stokes and anti-Stokes peaks. 
The vertical scale is in arbitrary units. 



V. DISCUSSION 

For simple systems, those (for example) with a single Ra- 
man transition driven by relatively long pulses, it will usu- 
ally be most efficient to continue using a multi-field model. 
Single-field simulations require very fine time-resolution, so 
they are computationally expensive for pulses with many op- 
tical cycles. The spectral range of the numerical field is cor- 
respondingly broad, typically covering many Stokes and anti- 
Stokes lines. 

In more complex situations, however, the single-field ap- 
proach will outperform its multi-field counterpart. For exam- 
ple, if a Raman interaction is probed by a beam that does not 



lie on the frequency comb defined by the pumping beams (e.g. 
as in [|l3|]), the multi-field approach will become much more 
complicated to implement. It will be necessary to define sepa- 
rate arrays for the pump and probe Raman "ladders" of Stokes 
and anti-Stokes lines, an issue that we avoided in section |^ by 
replacing the pump stage of the process with an initial condi- 
tion for the polarization. With a single-field model, the probe 
pulse and its Raman sidebands simply get superimposed on 
the overall spectrum, where they will be offset from the fre- 
quency ladder defined by the pump beams. 

Another situation in which the multi-field model will run 
into difficulty is where there are multiple Raman resonances. 
Although the treatment in this paper has been restricted to 
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FIG. 2: lOfs probe pulse incident on a medium with an initial polar- 
ization of pi2 = 0.025;. Here we compares three simulation results: 
(a) single-field GFEA simulation, (b) multi-field simulation, and (c) 
single-field SVEA simulation. The dashed lines help compare the 
relative heights of the first Stokes and anti-Stokes peaks. The verti- 
cal scale is in arbitrary units. 



a simple two-level Bloch equation description of the Raman 
medium, additional Bloch equations can easily be added, even 



if there are coupled multi-level interactions (as for example in 
ijlj]). It is only necessary to describe those transitions appro- 
priately, and to modify the polarization terms acting on the 
propagating field. This procedure is considerably more diffi- 
cult to handle in the multi-field case, which is based on field 
components separated by a particular Raman transition fre- 
quency. Additional Raman resonances complicate the theory; 
not only must extra detuning factors be added to the equa- 
tions, but it is also necessary to work out which field compo- 
nent is nearest to each new driving term. With a wideband 
single-field model, on the other hand, any new sidebands or 
resonance effects appear automatically in the spectrum, and 
no special measures need to be adopted to handle them. 

The usefulness of our single-field approach is not restricted 
to the Raman interaction described in this paper It is not just 
more easily extended to more complex Raman materials in- 
volving e.g. multiple transitions than the standard multi-field 
model. It would be equally valuable for a near-degenerate op- 
tical parametric oscillator, or indeed any system where two or 
more field components start to overlap as the pump or probe 
pulses get shorter. 



VI. CONCLUSION 

We have considered how best to model the multi-frequency 
field in wideband Raman generation experiments. Rather than 
using multiple field envelopes, with one at each Stokes or anti- 
Stokes frequency, we instead use a single wideband field en- 
velope. This requires that the field be propagated taking into 
account wideband effects, as described by either the SEWA 
theory of Brabec and Krausz [Q, or the more general GFEA 
of Kinsler and New . 

Our single-field approach has three crucial advantages. 
First, it includes more physics, even compared to a multi-field 
approach enhanced by adding GFEA corrections to the propa- 
gation of the field components. Secondly, it deals effortlessly 
with the complications of overlapping spectra that occur in the 
multi-field case. Thirdly, it allows for extra Raman transitions, 
and other molecular details to be included more easily than is 
possible for the multi-field model. 

All of these factors ensure that our wideband single-field 
model not only extends the regime in which envelope-based 
methods can be utilised; but is also more versatile and places 
fewer restrictions on the kinds of Raman media that can be 
easily described. 
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